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TECHNICAL NOTE 3893 

THEORY AND DESIGN OF A PNEUMATIC TEMPERATURE PROBE AND EXPERIMENTAL 
RESULTS OBTAINED IN A HIGH- TEMPERATURE GAS STREAM 
By Frederick S. Simmons and George E. Glawe 


SUMMARY 

The basic theory of pneumatic temperature probes and deviations 
from this basic theory in practical application are presented. Design 
requirements and operating considerations are discussed. Temperatures 
between 1600° and 4000° R were measured with an experimental probe in 
hydrocarbon combustion exhaust gases and compared with those measured by 
thermocouples and a line-reversal pyrometer. Reasonable agreement among 
the instruments was obtained. 


INTRODUCTION 

An accurate temperature measurement of a gas stream becomes increas- 
ingly difficult with increasing temperature, and most of the better known 
methods have some shortcomings. For example, immersion instruments such 
as thermocouples or resistance thermometers that are capable of with- 
standing high temperatures require large corrections of the indications 
to account for heat-transfer and recovery effects. Optical techniques, 
such as sodium- line reversal, usually measure some mean temperature along 
the optical path and require the presence of chemicals in the gas stream. 

A method which appears to offer several advantages is the use of a 
pneumatic temperature probe. This instrument consists essentially of two 
restrictions through which a sample of the hot gas is successively passed 
and between which the gas is cooled. From continuity of mass, the stream 
temperature is a function of the pressure measured at each restriction, 
the temperature at the downstream restriction, and the areas and discharge 
coefficients of the restrictions. Thus, a thermodynamic measurement of 
gas temperature can be obtained beyond the range of most commonly used 
thermocouple materials. 

Publications (through 1951) on pneumatic -probe methods of measuring 
gas temperature are listed in reference 1, and later work is described in 
references 2 and 3. This report reviews some of the considerations 
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entering into pneumatic -probe design and compares a pneumatic probe with 
thermocouple pyrometers and a line-reversal pyrometer. The pneumatic 
probe herein described used a nozzle for the upstream restriction and a 
venturi tube for the downstream restriction, with sonic flow through both 
restrictions. Measurements were made in the temperature range from 1600° 
to 4000° R, at subsonic stream Mach numbers, and at static pressures of 
the order of 1 atmosphere. This work is part of the research program in 
high-temperature measurements being conducted at the NACA Lewis laboratory. 


THEORY 

Basic Equations 

A schematic diagram of the pneumatic probe with pressure and tem- 
perature notation is shown in figure 1. (Symbols are defined in 
appendix A . ) 

The basic equation for mass-flow rate, which assumes that a perfect 
gas is expanded isentropically and one -dimensionally in an orifice, is 
given by 


w 


PA / 2 

a/t VR 



( 1 ) 


Defining a discharge coefficient C as the ratio of the actual 
mass-flow rate to the ideal mass-flow rate and equating the actual mass- 
flow rates of the two restrictions in series yield, for constant R: 




Solving for the temperature ahead of the front nozzle gives 
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where 



Application of equation ( 3 a) to a pneumatic temperature probe in 
general involves the measurement of the pressures Pq, P2> P3, and P4 
and the temperature T3; calculation of the quantities Fp and F3 from 
the measured pressures with knowledge of the effective specific-heat 
ratios Xq 2 an< ^ T3 4> an< ^ knowledge of the ratios Cg/C^ and Ag/A^. 


For the special case where the flow through both restrictions is 
sonic , it can be shown that 


where 



f sir) 



( 3 b) 


If it is assumed that C2A2/C4A4 is a constant and To 2 = V3 4; equa- 
tion ( 3 b) reduces to 


Tq = Constant * 



( 3 c) 


Effects of Heat Transfer on Nozzle Discharge Coefficient 

The quantity (C2A2/C4A4) 2 in equation ( 3 b) is usually considered 
as a constant. In practice, there are variations in this quantity 
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resulting from changes in the physical areas of the nozzles and variation 
in the discharge coefficients with Reynolds number and heat transfer 
between the gas and the nozzle surfaces. It is difficult to determine 
(A 2 /A 4) 2 accurately by direct physical measurement of the nozzle throats. 


since this quantity is proportional to the fourth power of the diameters. 
However, from equation (3a) it can be seen that 


CsA2\ 2 

, C 4 A 4/t 0=T3 


P3 F3\ 2 

P n Ft / 
v0 1 / Tq^T 3 


Therefore, the quantity (C 2 A 2 /C 4 A 4 ) 2 can be obtained by a calibra- 
tion, that is, operation of the probe at room temperature where Tq = T 3 . 
The factor (C 2 /C 4) 2 can then be computed from known relations between 
discharge coefficients and Reynolds numbers and separated from the term 
(C 2 A 2 /C 4 A 4 ) 2 . This development thus yields an effective area factor 
(A 2 /A 4) 2 which in application would be altered only by thermal expansion 

or erosion and particle deposition. The effect of thermal expansion for 
a hot run can be calculated with knowledge of nozzle temperatures or can 
be reduced to a negligible quantity by proper cooling of the restrictions. 
The effects of erosion or particle deposition can be accounted for by 
periodic calibration at room temperature. 

The quantity (C 2 /C4 ) 2 to be applied to the solution of Tq for a 

hot run can be treated by using the expression suggested in reference 4, 
which shows that the discharge coefficient of a flow nozzle, when the 
fluid is a gas and there is heat transfer between gas and nozzle, is 


C - 1 




+ 


1 

4 



to 


where t/t n is the ratio of gas to nozzle-wall temperature, Z the 
length of straight throat section, Z ' the axial length of the rounded 
approach section, the nozzle throat diameter, and 



The calculation of C 2 involves the Reynolds number of the upstream 
nozzle; this Reynolds number is a function of the temperature of the gas 
which is being measured. An iterative calculation would, therefore, be 


5040 


NACA TN 3893 


5 


necessary, t>ut an adequate value of Tq f>or "this Reynolds number can be 
obtained from equation (3c) using the value of (C 2 A 2 /C 4 A 4) 2 obtained in 
a cold calibration. The computations are facilitated by the use of fig- 
ure 2, which shows the Reynolds number for critical flow through a 
l/ 8 - inch- diameter nozzle as a function of gas temperature and pressure, 
and figure 3, which is a plot of equation (4) for the nozzle used 
(D = 21 = l') . 


Effective Specific-Heat Ratio 

Equation (l) was derived on the assumption that the ratio of spe- 
cific heats of the gas is a quantity independent of the process the gas 
undergoes. At room temperature this is true for real gases, but as the 
temperature is increased, energy absorbed by internal vibration of the 
molecules becomes increasingly important. Consequently, the specific 
heats and their ratio become functions of temperature. However, the 
transfer. of energy into the vibrational mode does not occur with the 
same rapidity as that among the rotational and translational modes, the 
reason being that only certain types of collisions, which occur relatively 
infrequently, can excite a molecule in vibration. A characteristic or 
relaxation time for the vibrational mode of the gas may be defined and 
shown to be a function of the pressure, temperature, and composition of 
the gas. Thus, should a gas undergo a change, for instance, an expansion, 
in a time much shorter than the vibrational relaxation time, the process 
would correspond to one in which y retains a value corresponding to the 
absence of any vibrational degrees of freedom. Reference 5 is concerned 
with this effect in flow processes involving hydrocarbon products. 

Three other effects further complicate the problem of the nozzle 
flow processes occurring in a pneumatic probe: 

• 

(1) If the pneumatic probe is immersed in a high-velocity gas 
stream, a compression and an expansion, not necessarily of like duration, 
occur at the upstream nozzle. 

(2) If the gas stream is of sufficiently high temperature that dis- 
sociation is appreciable, the cooling of the gas within the probe would 
involve an uncertain degree of recombination with resultant uncertainty 
in the specif ic -heat ratio and the gas constant. 

(3) Most gases of interest are mixtures; a net effect would involve 
combinations and interactions of the preceding effects among the 
components . 

An approximation for the case where the compression ahead of the 
upstream nozzle occurs in a time of the same order as the relaxation 
time for the hydrocarbon gases is presented in appendix B. A relation 
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is established wherein an effective y"q 2 -*- s re l a ted to equilibrium 

high-temperature To; equilibrium low-temperature T3 4^ and free -stream 
Mach number. 


PROBE DESCRIPTION AND DESIGN CONSIDERATIONS 
Probe Description 

A photograph of the experimental probe is shown in figure 4; details 
of the first-nozzle assembly, which had a l/8-inch-diameter throat, are 
shown in figure 5. The second nozzle was geometrically similar with the 
exception of a diffuser section, which in effect made it a venturi tube. 
Both nozzles were constructed from a chrome -copper alloy. The probe was 
designed for testing at gas-stream total pressures of the order of 1 
atmosphere to utilize the large-capacity laboratory exhaust system (ap- 
proximate pressure, 10 in. Hg abs) for pumping. A diameter ratio 
^th,2/®th,4 of unity was chosen to ensure sonic flow in both nozzles at 
the anticipated range of high-temperature operation. The thermocouple 
for measuring gas temperature T3 was chromel-alumel, 0.010 inch in diam- 
eter, and had an immersion length of 50 diameters. 


Design Considerations 

Choice of restriction shape . - The restriction configurations of 
the probe described were chosen because : 

(1) Some preliminary experiments in exhaust gas streams containing 
a relatively large amount of unburned solid particles had shown that 
throat area and discharge coefficient were changed less by erosion and by 
dirt deposition in the case of flow nozzles than in the case of sharp- 
edged orifices. 

(2) It was felt that constancy of throat area could best be main- 
tained if the restriction were cooled and that this cooling could be 
effected more easily for a thin-walled nozzle than for a sharp-edged 
orifice. 

(3) At constant Reynolds number, the discharge coefficient of the 
nozzle remains constant with pressure ratio in supercritical flow (ref. 

6), whereas that of the sharp-edged orifice does not. 

(4) Pressure loss at the downstream restriction was least for the 
venturi-type design; thus, severity of pumping requirements was minimized. 
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The disadvantages of the nozzle are: 

(1) The effect of heat transfer from the gas to the cooled surface 
of the nozzle is greater than that for the sharp-edged orifice. This 
heat transfer affects the boundary -layer flow and hence the discharge 
coefficient. 

(2) Mechanical construction of a nozzle is more difficult than that 
of a sharp-edged orifice. 

(3) For throat diameters in the range of 0.04 to 0.1 inch, the var- 
iation of discharge coefficient with Reynolds number at constant pressure 
ratio is greater for a nozzle than it is for a sharp-edged orifice. 

The advantages itemized for the nozzle were considered to outweigh 
the disadvantages, since an approximate theory was available to correct 
for the effect of heat transfer on the discharge coefficient (ref. 4 ) . 

The size of the throat diameter chosen was based upon a compromise 
between pumping requirements and the advantage of a large throat area in 
achieving high Reynolds numbers and minimizing effects of dirt deposition. 

Pumping requirements . - Pumping requirements may be estimated by 
using equation (3a) and assuming C 2 /C 4 = F^/F^ = 1 and by choosing a 
convenient temperature T 3 . This temperature is chosen high enough to 
be well above the dew point and low enough to be measured accurately. 

Then equation (3a) becomes 


Assuming sonic flow in the first nozzle and slight pressure loss 
between stations 2 and 3 gives 

£3 < £2 w 1 


so that the required diameter ratio is specified by the inequality 


Assuming sonic flow in the second nozzle and 60 -percent recovery of 
the velocity head in the diffuser, the approximate pressure ratio at the 
second nozzle is 



( 5 ) 



(?) 
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( 8 ) 


Therefore , 


P, 


5 


P 


3 0.8 


( 9 ) 
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If sonic flow exists also in the first nozzle, equations (6) and (9) 
yield 


This inequality establishes the pumping pressure that must he maintained 
at the downstream end of the probe to achieve sonic flow in both nozzles 
at high-temperature operation. Availability of adequate pumping facili- 
ties to maintain this over-all pressure ratio would eliminate the neces- 
sity of measuring static pressures pg and p^, which need be known only 
if sonic flow does not exist in both nozzles. From equation (7) it can 
be seen that when Tq = T 3 the diameter ratio 2/-9th 4 must be equal 

to or less than 0.7 for sonic flow to exist in both nozzles. However, 
for this case the pumping requirements would be more severe at high- 
temperature operation. For a diameter ratio of unity, subsonic flow will 
exist in the first nozzle for a calibration at Tq = T3, but sonic flow 
in both nozzles at high -temperature operation can then be attained with 
less severe pumping capacity. 

Cooling requirements . - Conventional heat -transfer formulas are 
adequate for probe cooling calculations. By arbitrarily selecting some 
fixed dimensions and from knowledge of anticipated operating conditions, 
the dependent quantities such as coolant flow rate and pressure drop can 
be found. A detailed procedure for calculation is presented in appendix 
C. 


The use of the pneumatic probe involves, in the most general case, 
the measurement of four pressures, one gas temperature, and two nozzle 
wall temperatures. 

The total pressure Pq of the gas stream may be measured by any of 
three approaches: A separate total -pres sure probe can be used; the flow 

in the probe can be stopped, and the P3 gage used to make the measure- 
ment; or a pressure tap can be provided on the face of the upstream 
nozzle (P^) . In this last method, however, the total pressure P^ so 



( 10 ) 


Pressure and Temperature Measurements 
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measured was not fully equal to Pq, and varied primarily with Mach number, 
so that a pressure calibration relating P-^/Pq to Mach number M was 

required (fig. 6 ). The measurement of P 0 with a separate probe , when- 
ever possible, appears to be the most advantageous method. Measurement 
by the P]_ tap not only has the defect previously noted, but also suffers 
from the fact that small tubing must be used which penalizes the response 
time. The method of stopping the flow in the probe has the disadvantage 
that a gage would have to be used which has a greater range than 

would be necessary for the measurement of under flow conditions. 

Also, the possibility of internal condensation would present the problem 
of time necessary to reestablish continuity of mass when the pumping is 
resumed. 

Static pressures p^ and p^ need be measured only if there is any 
uncertainty that sonic flow exists in both nozzles. Measurements made 
with the p 2 and p^ taps on the present probe confirmed that equations 
(7) and (10) were reasonable approximations and sonic flow existed at 
both restrictions. 

The measurement of total pressure P 3 upstream of the second nozzle 
can be obtained with a wall tap if the duct area upstream of the second 
nozzle is sufficiently large to substantially constitute a stagnation 
region. 

The region upstream of the second nozzle is then also a convenient 
location for measurement of total temperature T 3 by means of a thermo- 
couple. This temperature should also be sufficiently low so that the 
radiation and conduction errors of the thermocouple are negligible. 

Knowledge of the nozzle wall temperature was necessary for the com- 
putation of discharge coefficient under the condition of heat transfer 
(ref. 4). Experiments with some preliminary probe designs had shown 
that, with the cooling arrangement shown in figure 5, the nozzle wall 
temperature did not greatly exceed the water temperature. Hence, direct 
measurement of nozzle wall temperature by thermocouples was not necessary. 


APPARATUS 

High-Temperature Tunnel 

The tests were performed in the high-temperature water-cooled tunnel 
described in reference 7. This facility is capable of continuous opera- 
tion at temperatures from 1400° to 3700° R, Mach numbers from 0.2 to 
1.0, and pressures from l/2 to 2 atmospheres. Using oxygen, operation 
for short periods to a temperature of 4000° R is possible within re- 
stricted ranges of velocity and pressure. The fuel throughout these 
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tests was 60-octane gasolene. Total pressure was measured with a water- 
cooled probe upstream of the 4- inch nozzle, and static pressure by a meas- 
urement of receiver pressure. 


Comparison Instruments 

An aspirated chromel-alumel thermocouple probe was used for compari- 
son from 1600° to 2500° R. Reference 8 shows data indicating this probe 
to be sensibly free from radiation and conduction errors in the range of 
conditions encountered in these tests. 

Two bare -wire, butt-welded, platinum- 13-percent -rhodium - platinum 
thermocouples with water-cooled supports were also used for comparison. 

One platinum thermocouple was a 0.020-inch-diameter wire with an immersion 
length of 50 diameters, and the other platinum thermocouple was a 0.032- 
inch-diameter wire with an immersion length of 28 diameters. The junc- 
tions of these thermocouples were butt -welded and then swaged to reduce 
the welding "bead" to wire diameter so that the configuration represented 
a cylinder in crossflow with uniform diameter. The sodium-line -reversal 
instrument was essentially a self-balancing adaptation of the standard 
instrument and is described in detail in reference 9. 


TESTS ANR RESULTS 
Calibration Tests 

The calibration coefficient (C 2 A 2 /C 4 A 4) 2 was experimentally ob- 
tained using room-temperature air, at total pressures from 20 to 50 
inches of mercury absolute, and at Mach numbers from 0 to 1.0. The value 
obtained was constant to within the experimental accuracy. The mean 
value of (Ag/A ^.) 2 was found to be 0.842. The nozzle -throat diameters 
were specified as being equal, so theoretically this factor should be 
unity. The calibration therefore indicates that the actual nozzle diam- 
eters were approximately within 2 percent of specified diameter assuming 
one oversize and the other undersize. These results indicate the impor- 
tance of a calibration, since an error of about 16 percent in Tq would 
result from assuming (Ag/A ^) 2 to be unity. 


High-Temperature Tests and Results 

The experimental pneumatic probe was tested from 1600° to 4000° R 
at a stream total pressure of 1 atmosphere and at Mach numbers of the 
order of 0.5 in overlapping ranges of the comparison instruments previ- 
ously described. 
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To eliminate effects of gradients in temperature and pressure, the 
comparison "between probes was made by moving each probe, in turn, to the 
same location in the gas stream with tunnel conditions held constant. 


For the line-reversal measurement, sodium bicarbonate powder was 
injected into the gas stream through a water-cooled probe upstream of the 
test region. This method of injection provided a local stream of sodium 
vapor in the test section so that a localized temperature was obtained. 
During injection of the compound, the other instruments were retracted 
to avoid particle deposition. 

During the tests, the temperature ahead of the second nozzle T3 
was of the order of 700° R. The pneumatic -probe temperatures were cal- 
culated by equation (3b) with the following assumptions: 

(1) The effective specific-heat ratio Yq 2 can calculated by 
the method presented in appendix B. 

(2) The effect of heat transfer on the nozzle discharge coefficients 
can be calculated by equation (4) with nozzle wall temperature taken as 
the mean cooling-water temperature. 


The radiation, conduction, and recovery corrections of the platinum 
thermocouple probes were calculated by using methods and values given in 
references 7 and 10. The maximum total correction for the 0.020-inch- 
diameter platinum thermocouple was 270° F, and the maximum total correc- 
tion for the 0.032- inch-diameter platinum thermocouple was 360° F. The 
results of these comparisons are plotted in figure 7(a). At 3700° R the 
agreement with comparison instruments was of the order of ±2 percent. 

These data for the pneumatic probe were also computed by using the 
calibration value of (C2A2/C4A4.)^ in equation (3b), in order to show 
the result of neglecting the heat-transfer and Reynolds number effect on 
the coefficient of discharge of the nozzles. The results of this com- 
putation are in figure 7(b), and they show a systematic error. 


In tests made at Mach numbers from 0.3 to 1.0 at a stream total 
temperature of 2460° R the comparison instrument was the aspirated chromel- 
alumel thermocouple. The results of these tests are shown in figure 8. 

The squares represent temperatures calculated using T"q 2 = l"i 2 = ^ 3 4 
and show a marked dependence on stream Mach number. From this result 
equation (B5) was deduced. Temperatures calculated using equations (B5) 
and (B6) are indicated in the same figure by the circles. It is noted 
that the corrections bring the pneumatic -probe data to within about 1 
percent of the comparison instrument. A sample calculation is given in 
appendix E. 
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The gas stream was sufficiently clean that systematic changes in 

the coefficient (C 2 A 2 /C 4 A 4 ) 2 as determined hy calibrations with room- 

temperature air before and after each series of tests were less than 1 
percent. The average duration of exposure was about 1 hour. 


DISCUSSION 

In view of the lack of a means of accurately establishing the true 
gas temperature in these tests, the absolute accuracy, as such, of the 
pneumatic probe cannot be stated. However, the agreement among the 
methods, as shown in figure 7(a), is fairly satisfactory from 1600° to 
3800° R. 

Above 3800° R the pneumatic probe indicates considerably higher 
temperatures than those of the sodium-D-line-reversal instrument. A 
possible cause is the fact that these data were taken using oxygen in 
addition to air in the combustion. The pneumat ic -probe temperatures were 
calculated using the same assumptions as for hydrocarbon-air combustion. 
Not only are these assumptions questionable, but at these temperatures 
the effects of dissociated water vapor and carbon dioxide in the products 
become apparent. A detailed analysis of these factors was considered to 
be beyond the scope of the present work. Furthermore, from the experi- 
mental procedure, it is felt that errors exist in the D-line measurements 
at these higher temperatures such as those due to sodium contamination 
(ref. 9) . The experimental procedure involved the addition of oxygen to 
a very rich fuel-air mixture while observing the indication of a platinum 
thermocouple located in the cooler outer boundary of the jet. While the 
thermocouple indication and the pneumatic -probe temperature increased 
uniformly with increasing amount of oxygen, the D-line temperature showed 
little or no increase. 

The data shown in figure 8 indicate that errors of the order of 5 
percent may result from not considering the compression ahead of the up- 
stream nozzle when the probe is immersed in a gas stream of Mach number 
near unity. The empirical relations represented by equations (B5) and 
(B6) were selected from these data, which were obtained at a fairly 
accurately measurable gas temperature of 2460° R. 

The use of equations (4) and (B6) which take into account heat trans- 
fer and effects of flow variations, respectively, in the nozzles may be 
considered as corrections to the basic equation (3c) of the pneumatic 
probe. Because the corrections are, at most, a few percent in magnitude, 
the validity of equations (4) and (B6) is not clearly shown by the ex- 
perimental results, particularly in the instances where the corrections 
tend to counteract one another, such being the case in the sample calcu- 
lation. On the other hand, since the corrections are small, the use of 
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crude assumptions in the derivation of the equations is justified in 
that large errors may he present in the corrections without introducing 
significant errors in the temperature measurements. 

Appendix D shows orders of magnitude of the various possible errors 
mentioned throughout the text. 


CONCLUDING REMARKS 

An experimental probe based on the design criteria herein reported 
was operated to 4000° R and compared with other types of pyrometers. 
Reasonable agreement among the instruments was obtained. The following 
recommendations based on this investigation are made: 

1. Flow nozzles should be chosen in preference to thin-plate orifices. 

2. Temperature ahead of the downstream nozzle should be low enough 
to be accurately measured, but above the dewpoint of the gas at the 
pressure ahead of the downstream nozzle. 

3. Sonic velocity should be maintained in both nozzles when the 
probe is used to obtain free-stream total temperature. A nozzle diameter 
ratio of unity can be used to decrease pumping requirements, in which 
case the room-temperature calibration must accept subsonic flow in the 
first nozzle. 

4. Free-stream total pressure should be measured with a separate 
total-pressure probe if possible. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 12, 1956 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 
cross-sectional area of nozzle throat 
discharge coefficient 

specific heat of gas at constant pressure 

specific heat of water 

diameter 



function 

acceleration due to gravity 
thermal conductivity 
probe length 

length of straight section of nozzle 

axial length of approach section of nozzle 

Mach number 

total pressure 

static pressure 

heat- transfer rate 
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R specific gas constant 
Re Reynolds number 
T total temperature 

t static temperature 

W fraction of water vapor in gas 

w ideal mass -flow rate 

x wall thickness 

y radial depth 

X specific-heat ratio 

p density of gas 

\i viscosity of gas 

t relaxation time of gas 

Subscripts: 
e external 

g gas 

i internal 

n wall at nozzle throat 

r reference 

s sonic 

th throat 

w cooling water 

0 free stream 

1 immediately ahead of upstream nozzle 

2 throat of upstream nozzle 
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3 ahead of downstream nozzle 

4 throat of downstream nozzle 

5 behind downstream nozzle 
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APPENDIX B 


- 


ESTIMATION OF EFFECTIVE SPECIFIC- HEAT RATIO 

Effective values of the specific-heat ratio for the flow process at 
the upstream nozzle may be defined by the following equations which assume 
the velocity at station (l) to be zero and the flow to be one- dimensional: 


T r 


= 1 + 


*0 


4 




T i r 0 ,i - 1 z 


^i Ti 2 ~ 1 ? 


(Bla) 


J 


The temperature ratio for the entire process is described by 


T- = 1 + 
t 2 


ro, 2 - 1 


(Bib) 


where the effective specific-heat ratio of the entire process Xq 2 
obtained by simultaneous solution of equations (Bl) 


r °’ z ' 1 + 5| 


i + 


fl,2 " 1 


MgMl + 


■^O 


4 




1 + ''V 1 

which, for < 1 can be approximated by 


- 1 


r o,2 ■ U,2 - |(r 0 ,i - r 0 ) ^ 

For sonic flow in the first nozzle (Mg = l), 

r 0 ,2 - ?i,2 - (r 0 ,i - r 0 ) 4 


(B2) 


(B3) 


18 


NACA TN 3893 


Equation (B3) has significance only for the case where the expansion 
through the nozzle occurs rapidly compared with the relaxation time of 
the gas, in which case Tq 2 * T 3 4 > and - To 1 ma y range from Xq to 

Y± 2’ ^Pending on the stream Mach number Mq and the probe geometry. 

For hydrocarbon combustion gases where nitrogen is the main con- 
stituent, reference 11 shows the relaxation time to be 

x oc Vt/( P W) (B4a) 


or 


X 



*r 

w 


(B4b) 


where the subscript r indicates some reference conditions. Reference 

11 gives data showing x = 2.3xl0 “ 4 second for t = 1034° R, p = 1 
atmosphere, and W = 0.02 for a nitrogen - water- vapor mixture, so that 

T * 1 . 4x10” 7 Vt/(pW) (B4c) 

where p is in atmospheres. 

If the velocity changes in the compression and expansion occurring 
at the upstream nozzle are assumed to be linear with and to occur within 
distances of the order of the probe diameter D e , the order of magnitude 

of gas velocity at which relaxation effects become significant is 2D e /x # 

For a probe diameter of l/2 inch, this is roughly equivalent to stream 
Mach numbers of 0.2 and 0.4 for temperatures of 1000° and 4000° R and a 
pressure of 1 atmosphere for hydrocarbon combustion products. Thus, for 
this probe size with M 2 = 1 , n 2 may be ^ aken as the room- temperature 

specific-heat ratio for the particular gas mixture) it is a function of 
the composition of the gases, which, of course, varies with the tempera- 
ture attained in the combustion. The specific-heat ratio To 1 is a 

function of the equilibrium specific-heat ratio y*q, which varies with 
temperature and composition, and of the stream Mach number Mq. 

On the basis of test results at a constant temperature with a varia- 
tion in Mach number (fig. 8 ), the following empirical relation for To 1 
was established: 


* 0,1 = r 0 + (tt 3 ,4 - r 0 ) Mo 


(B5) 
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Substitution of this relation in equation (B3) yields the following 
evaluation for r 0 , 2 ln terms of y Q , r 3 ^ 4 , and Mq: 


^0,2 - r 3 ,4 - (r 3 , 4 - r 0 ) 


(B6) 


Application of equation (B6) in conjunction with equation (3) may 
be facilitated by use of graphs such as figures 9, 10, and 11, where r 0 

and T3 4 are given as functions of gas temperature. o as a function 
y r / no 

f s (t 0 2 ) 2 

of* Xq and Mq, and y as a function of Tq^ and 

Again an estimation of the temperature of the gas stream should be 
adequate. 
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APPENDIX C 

COOLING REQUIREMENTS 

The water- cooling requirements for the probe can be estimated from 
four heat- transfer-rate relations and one pressure- loss relation. 
Nomenclature used is illustrated in figure 12. 

(1) The rate of heat transfer from the gas stream to the outer shell 
of the probe q e can be obtained from a specification of probe diameter 

D e , exposed length L e , gas pressure, temperature, and velocity, and 
outer- shell temperature T e . For this computation, gas temperature may 
be taken as Tq and outer- shell temperature as average cooling- water 
temperature T w av . It is specified that the cooling- water temperature 
always be well below the boiling point, say 160° F. Then 

q e = 3tD e L e h e (TQ - T w ^ av ) (Cl) 

where h e may be obtained from 

^-0.24 (Re D )°- 6 (02) 

K g e 

(ref. 12, p. 222). 

(2) The rate of heat transfer from the aspirated gas to the inner 
shell of the probe q^ can be obtained from the specification that 

q i = c P,g( T 0 " T 3) w g ( C3 ) 

This rate can be related to the exposed internal length L^, internal 
diameter D^, Reynolds number Re^, and inner- shell temperature T^ by 

T - T 

!i = c p,g(^0 “ ^3) v g = ki^i^i /Tq - T • \ (C4) 

in w^) 

where h^ may be obtained from 

HT 1 = 0.02(Re D ) 0 ’ 8 (C5) 

*g 1 

(ref. 12, p. 170). For this computation, may also be taken as equal 
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In equations (C2) and (C5), the gas viscosity and thermal conductivity 
may be evaluated at any convenient fixed gas temperature such as 
(T 0 + T3 ) /2 , since the heat- transfer coefficient varies only as the 0.2 
power of the gas temperature. 

(3) The total rate of heat transfer to the water is the sum of q 
and q^ and produces a temperature rise AT W . 

q e + 9i = «w c w^w (C6) 

For an assumed inlet water temperature of 60° F and a specified maximum 
water temperature of 160° F, 

AT W » 100° F 

(4) The pressure difference needed to produce the required water flow 
rate will be computed under the arbitrary simplifying assumptions that a 
linear pressure drop exists in the cooling annulus, that pressure loss at 
the 180° reversal of flow direction at the probe tip can be neglected, 

and that the annulus may be considered as having a mean diameter 
D av = ( D e + D i)/ 2 J a length 21^, and a radial depth of y = l/2 

[(D e - D. jJ/2 - 3x], where x is the average wall thickness of the three 
concentric tubes shown in figure 12. Then the applicable equation iss 


A P W 


24 ^w L i 

*Pwy 3D av ” w 


(C7) 


(ref. 12, p. 124). 

If D e , L g , L^, x, Tq, T3, T w ^ ay , AT W , and Wg are chosen in ad- 
vance, then equations (Cl) and (C3) give q e and q ± directly for use 
in equations (C6), equations (C4) and (C5) yield the required D^, equation 
(c6) yields the required w w , and equation (C7) yields the approximate 
A Pw* 

This cooling must be accomplished within the following limitations: 

(1) Water pressure available 

(2) Heat-transfer rate obtainable between the water and the outer 
and inner probe walls 

(3) Heat- transfer rate obtainable between the inner probe wall and 
the aspirated gas 
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( 4 ) Local water-passage temperatures low enough not to allow forma- 
tion of steam pockets 

In the event that some of the requirements as outlined are not compatible 
with a probe design and application, a measured temperature T 3 may he 

inconsistent with the recommendation that it he at an accurately measura- 
ble temperature and above the dew point. In this event a heat exchanger 
can be placed upstream of station 3 to bring T 3 within a proper range. 
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APPENDIX D 


ORDERS OF MAGNITUDE OF POSSIBLE ERRORS 
Effect of Heat Transfer and Variation of Reynolds 
Number Using Flow Nozzles 

In the sample calculation given in appendix E, it is seen that neglect 
of the effect of heat transfer and the variation of Reynolds number be- 
tween calibration with room-temperature air and application at high tem- 
perature would have resulted in an error of -2^ percent. In the present 

tests with hydrocarbon gases from 1500° to 4000° R, such potential errors 
were less than 4 percent and greater than 2 percent. In the sample cal- 
culation, had the Reynolds number variation only been considered, from 
figure 3, it can be seen that an error of +5 percent would have resulted. 
These conclusions, of course, assume the validity of equation (4). 


Uncertainty in Effective Specif ic -Heat Ratio 

From figure 11 it can be seen that the error in temperature measure- 
ment when a value of 1.4 is assumed for the effective specif ic -heat ratio 
I"0,2 ranges from 2.6 to 7.8 percent for actual values of Pq g from 

1.35 to 1.25. The actual values of p Q , 2 > as discussed previously, depend 

in a complex manner on the gas properties, composition, and velocity, and 
on nozzle geometry and method of total-pressure measurement. 

In the present tests, the largest error due to neglect of the effect 
of stream velocity shown in figure 8 would be 4.4 percent, the datum 
obtained at M a 1, Tq * 2460° R. It would not be unreasonable to con- 
clude, however, that for hydrocarbon-air combustion gas streams of total 
pressure of the order of 1 atmosphere and Mach number less than 0.5, 
negligible errors result from the assumption that r - y 

'0,2 '3.4* 


Uncertainties in Internal Temperature and Pressure Measurements 

As shown by equation (3), errors in the measurement of T 3 are 

directly reflected in the calculation of T Q . Such errors are related 

to the temperature, pressure, and gas velocity at station 3 and to the 
diameter and immersion length of the thermocouple at that station. With 
reasonable attention to the design of the probe and the thermocouple 
such errors may be shown to be negligible (refs. 10 and 13). 
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Errors in the measurements of pressures P Q and may he con- 

sidered negligible, and if pumping facilities are adequate to produce 
sonic flow in both nozzles for high-temperature operation, it is not 
necessary to measure p 2 and p 4 . However, in calibration at room tem- 
perature, a probe of unity nozzle -diameter ratio will have subsonic flow 
in the front nozzle; therefore, static pressure p 2 must be accurately 

known. Hence, it is important that the static tap be properly located, 
free of burrs, and so forth. A l/4-percent error in p 2 will reflect 

as a l/2 -percent error in Tq. 


Sampling Errors in Fluctuating Flow 

Reference 14 shows that errors in the mean temperature measurement 
for a rectangular wave-form fluctuating temperature may be calculated 
from 



where 1q n< ^ and T ac ^. are the indicated and actual temperatures, and 
and l-g are the fractional periods at which the gas is at tempera- 
tures a.nd t b- 

An example is given: For T^/^A = 4 8nd ^b/^A ~ 2 > T ind would 

be 12.5 percent higher than T act . Although this example is an extreme 

case, it is apparent that such errors in fluctuating flow measurements 
may be considerably higher than errors arising from the various other 
causes . 
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SAMPLE CALCULATION FOR EXPERIMENTAL PROBE 

Probe Calibration at Room Temperature 

mion .J+T prot> ® calibration at room temperature the following measured 
quantities and nozzle pressure ratios are specified: 

P l> p 0' PO' in * Hg aks pq IP 

p 2 , in. Hg abs * 

p 3- in - 24 .'4 0 

P 4 , in. Hg abs 

T 0' T 3' OR 

0.722 

P i^ P 3 0.501 

Hence, flow is subsonic in the upstream nozzle and sonic in the down- 
stream nozzle. Using equation (3a) with F-]_ = 0.442 and F 3 = 0.484 
gives 

M /P3 p 3\ 2 

w - \ml _ T - °- 842 

i 0 _1 3 

where 

M 2 

\ C 4/t 0 .T 3 

since Re 2 « Re 4 


Calculation of To for High-Temperature Operation 
with Sonic Flow in Both Nozzles 

In the calculation of Tq for high-temperature operation with 

sonic flow in both nozzles the following primary measured quantities 
are given: 
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p 0 , in. Hg abs 34.93 

P 3 , in. Hg abs 16.23 

T 3 > °R 


642 


The secondary measured quantities are 


p 0 , in. Hg abs 

T w> ^2,n> ^4,n> ^ 520 


Substitution of the given values into equation (3c) yields: 

T q - constant x t 3 (^j (0.842) (642) (yfj|§) 2 « 2500° R 

Correction for variation in discharge coefficient with Reynolds 
number and heat transfer . - From figure 2. Ren Lh . = 7.6xin3 ^ 

RS I>th,4 = 1 ' 8><104 * From 9, X Q = 1.286 and y^ „ = 1.374. The 


gas- to nozzle -temperature ratio is 

t2 T 0 


3,4 


u 2,n 

t 4 


1 + 


4) %n 


* 4.2 


T, 


4 ’ n 


1.0 


From figure 3, C 2 = 0.966 and C 4 = 0.955. Therefore, 


' _ /0.966 N 

^ C 4/Aiot v • 955 . 


= 1 . 022 . 


Correction for effective specific -heat ratio y . - Free-stream 

U y Cl 

pressure ratio is Pq/Pq = |^|| - 0.582. Therefore, Mq = 0.95. From 
figure 10, y* ^ 55 1.300; from figure 9, y = 1.374; from figure 11, 

r- , } 


F s(r 0 _ j 2) 


F s^3,4)j 


= 0.961. 
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Corrected Tq (eq. (3b)). - 


[M 

T & 

M 

F s(To,2) 

[_\v 

3 \ 3/ J 

\ C J 

F (r ) 

Lfs W 3,4 ; J 


T q = (2500) (1.022) (0.961) = 2460° R 
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Figure 1. - Schematic diagram of pneumatic probe. 
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Reynolds number per inch of mercury total 
pressure , Re/P , (in. Hg)' 1 
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Figure 2. - Reynolds number per inch of mercury total 
pressure as function of total temperature for critical 
flow in l/8- inch-diameter nozzle. Viscosity evaluated 
at total temperature; ratio of specific heats y, 1.35. 


5040 


Discharge coefficient 


5040 



MCA TN 3893 


CO 



First nozzle 


Figure 4. - Experimental pneumatic probe. 
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Figure 5. - Details of first-nozzle assembly. 
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Total -pressure error 
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Free-stream Mach number, Mq 

Figure 6. - Ratio of stream total pressure to indicated 
total pressure. Data obtained at pressures from 25 to 35 
inches of mercury absolute and temperatures from 2400° to 
3900° R. 
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(a) With consideration of heat-transfer effect on discharge coefficient. 


Figure 7. - Comparison of pneumatic -probe total temperatures with corrected indi- 
cations of other instruments. 
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(b) Neglecting heat-transfer and Reynolds number effects on discharge 
coefficient. 


Figure 7. - Concluded. Comparison of pneumatic -probe total temperatures 
with corrected indications of other instruments. 
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Figure 8. - Difference between calculated pneumatic- 
probe temperatures and aspirated-thermocouple tem- 
peratures as function of stream Mach number. Gas 
temperature, 2460° R. 
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Figure 9. - Equilibrium and corresponding low-temperature specific- 
heat ratios for high-temperature water-cooled tunnel. Low- 
temperature specific-heat ratios calculated at 700° R. 
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Figure 10. - Effective specific -heat ratio as function of equi- 
librium specific -heat ratio and stream Mach number. 
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Figure 11. - Variation of function 
TO, 2 and ^3,4- 


F s (TO, 2) 

F s(r 3 ' 4 ) 


with 
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Figure 12. - Schematic diagram for cooling calculations. 
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